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Abstract 

The simplest (3+l)D Regge calculus model (with three-dimensional discrete 
space and continuous time) is considered which describes evolution of the simplest 
closed two- tetrahedron piecewise flat manifold in the continuous time. The measure 
in the path integral which describes canonical quantisation of the model in terms 
of area bivectors and connections as independent variables is found. It is shown 
that selfdual-antiselfdual splitting of the variables simplifies the integral although 
does not admit complete separation of (anti-) self dual sector. 
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Path integral is probably the most efficient way to quantisation of the discrete theory. 
In particular, this approach is widely used in Regge calculus, including the physical 4D 
case [|l], 0. The path integral measure is usually chosen as the simplest among the 
invariant ones. Another approach is based on constructing state sum invariants of the 
4D Regge manifolds || |], [|] in analogy with how it is done in the 3D case where such 
invariants being functionals of only 2D triangulated boundary of such manifold exist 
and do not depend on the refinement of the interior tetrahedrons || 0, |§, || |10|. Thus, 
quantum measure is defined there by topological requirements. 

If description of the system is known in the canonical form, the most natural choice 
of path integral measure is that based on the canonical quantisation. Knowing canonical 
form of the theory requires that action for the system in the continuous time be known. 
In the given paper we consider path integral just for the simplest 3D Regge manifold 
evolving in the continuous time. Canonical form of this model as limiting case of the 



completely discrete 4D Regge manifold has been considered in [llfl in the area bivector 



and connection variables on the base of the canonical formulation of the general Regge 



manifold |L2] . The 3D Regge manifold consists of two similar tetrahedrons with mutually 
identified vertices. Remarkable feature of the related continuous time Regge calculus 
action is it's notational analogy and correspondence of constraints with the case of the 
continuum general relativity (GR) action in the Hilbert-Palatini (HP) form, whereas the 
number of the degrees of freedom is finite, and "large distances" (i.e. those considerably 
larger than the size of tetrahedron) are absent. Therefore the path integral for this 
system can be considered as that modelling path integral for the continuum GR on 
condition that long-distance effects (gravity waves) are excluded. The typical feature of 
the path integral measure obtained is that it vanishes for symmetrical configurations of 
the system. This indicates that quantum fluctuations at these points are definitive. 
Also path integral is convenient to handle the problem of implementing the so-called 



Ashtekar variables [13] in the Regge calculus (see review paper [14|) on quantum level. 
From general viewpoint, there are the two possible approaches to this problem: to put 
Ashtekar theory on the Regge lattice or to formulate Regge calculus in the Ashtekar-like 



variables. The recent approaches [Tj] seem to be of the first type, we try in our work 
the second one via splitting the antisymmetric tensor variables into the selfdual and 
antiselfdual parts. It turns out that although there is no complete separation of selfdual 
and antiselfdual sectors as it occurs in the continuum GR, the considered procedure 
simplifies path integral considerably. 

We start with the canonical structure of our model [IT]. The Lagrangian with all the 
constraints added reads: 

+ ^^(TT^OTT^+^^Tr^^+E^^*^- C 1 ) 

Here /i, u, A, p,... = 0, 1, 2, 3 run over four vertices of tetrahedrons; = —n£ a are 
bivectors of the triangles opposite (in 3D leaf) to the vertices /i; a, b, c, d,... = 0, 1, 
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2, 3 are local frame SO (3,1) indices (SO (4) in the case of Euclidean signature); f2" 6 C 
SO(3,l); AoB = (l/2)trAB for the bivectors A ab and B ab ; *A ab = (l/2)e abcd A cd ; A* B 
= A o (*B); h ab , h ab are multipliers at Gaussian constraints; the scalars enter also A 
= X)« U/xA^ and are, therefore, multipliers at the Hamiltonian constraints which are not 
Hp but Hn = Hp - A M Y,u>\ VVa^ ° tt A and 

A M 7r " = E £ M"W> (2) 
A,p 

H, = EMarcsin^p^ + ^X oA X ; (3) 

= E^A P (<^ + i[7T A ,^]), (4) 



X.p 

1 



-^V - 2 ^ (^mAa) (5) 

= n^, (6) 
it; = exp^x+VrO- (7) 

The completely antisymmetric symbol is defined so that e i23 = +1- The u v , (p", *tp" are 
nondynamical independent variables (not appearing in the kinetic term). These enter 
Lagrangian nonlinearly. The ip^v, and k^ v are variables entering L linearly of which 
ipfiv are dynamical ones. Of the six variables -?/v ( or OT ^y ^ wo are independent ones, 
and in JTT| the ip\2, ip23 (and k 12 , k 23 ) were chosen as the only nonzero ones. Now, in 



order to formulate the theory in the form explicitly symmetrical w.r.t. the permutation 
of the vertices 0, 1, 2, 3 we assume that all ipfiu, Kfiv are not identically zero, but are 
subject to the following four relations on ip^ u : 

^ = E = 0, (8) 

and four ones on n^ v : 

E = 0. (9) 

Let us take into account these relations on with the help of Lagrange multipliers 
Cfj,. Then we can consider 

^liv SIX independent variables, coefficients at the following 

constraints: 

EE TT^ * ^ + E CA = 0. (10) 

Xj^fi,u 

The are considered then as independent variables. 

Now passing to the Hamiltonian formalism we need the extended phase space co- 
ordinatised by the following canonical pairs (p,q): (ir^, fi^), (ip^ ((7^,(7^). The 
full system of constraints <3> consists of the first class rj and second class O ones. The 
coefficients at r\ in L are h, h and k^. The O consists of TC^, Z^ IU and the following 
constraints on the newly defined momenta: 

(> = (11) 
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and 



(12) 



Next we introduce the gauge conditions x = by the number of the first class 
constraints rj such that Det{i], x} 7^ 0. Then the path integral measure takes the standard 
form: 

dn x = e'S Ldt 8(® x )Det 1/2 {® x ,® x }DpDq 

= e i J Ldt Det 1/2 {e,e}5{Q)5{r ] )5(x)Vet{v,x}DpD(l (13) 

where <3> x = (<&, x)- For the specific form of the kinetic term for (p,q) = (ir, O) one 
gets DpDq = d 6 irV 6 tt, V 6 tt being the Haar measure on SO(3,l), and under the Poisson 
brackets (PB) one should imply the expression 

{/> 9} = n [U, flw] + U Mgn - ° n/n (14) 

(for the functions of purely 7r, fi; taking into account dependence on other variables is 
straightforward). In |ll], |12| this form was used as that describing Hamiltonian dynamics. 
Although quite expectable, this can be proved by considering standard kinetic term 
(l/2)trPf2 in the further extended phase space of canonical pairs (P, fl), (ip,ip), (C,C) 
and substituting into ([T3| ) instead of $ x the set $ x = ($ x ,t) extended by inclusion the 
constraints r specifying the 4x4 matrices P, Q to have the form Q e SO (3,1), P = 
fi7r, 7f = — 7r. Then Det{$ x , $ x } = Det{$ x , $ x }d( t ) where {•, -}d(t) are Dirac brackets 
(DB) w.r.t. the set r. The {•, -}d(t) turn out to coincide just with the earlier defined 
brackets (0). The Haar measure just arises when integrating over DPDQ the newly 
arising 5-functions of r, 

6(t) = n - 1)5^^ + P^). (15) 

A 1 

Now consider an important consequence of finiteness (in the Euclidean case) of the 
volume of the gauge group. The gauge transformations are some rotations changing 
Qfj,, *ip%. At the same time, the gauge fixing factor <5(x)Det{?7, x} m the path integral 
measure just arises upon dividing by the volume of the group of gauge transformations 
(generated by I class constraints rj via PB {rj, •}). Therefore, if we consider Euclidean 
case (or Lorentzian one in the sense of analytical continuation from the Euclidean case) 
we can insert this volume back into measure simply omitting gauge fixing factor without 
a risk of getting a new infinity. The measure takes the form 

dfi = e l I Ldt Bet 1/2 {Q, Q}5{Q)5{r])DpDq. (16) 

It is convenient to divide the set of second class constraints into two second class 
subsets $ and 1 so that Det{6, 6} = Det{$, $}r>u)- Let the set d consists of the four 
Hamiltonian constraints TC^, the 1 be the set of other second class constraints \& M , Z^, 
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C M , ifr^" To calculate the DB TC v }t>^) let us first define the projection TLj; of 
onto the subspace orthogonal w.r.t. the PB to the space of differentials of the set i: 

K = n »-H A ? z *> - E - E - E F ^ Xv (it) 

where the coefficients A, B, E, F are defined so that {7i^, l} = 0. For example, com- 
muting with constraints Z Xu and ^\ gives, respectively, 



- E B w + E ^ p ° ^ ^ * ^ + 7rA * = ( 18 ) 

and 

E = o. (19) 

These eqs. allow to find B^ v and Commutation with C v and \& Aiy allows to find 

A^ and i?^ from the analogous system. In the calculations of such kind commutators 
between different 7r-bilinears, 7r M o -n v and 7r M * tc u , arise. These commutators by ( |14|) 
reduce to the following trilinears: 

7T^ O [7T^, 7T A ] = 7T" * [tt", *7T A ] = 0, 

^*K,tt a ] = -V 2 5> PM „ A , (20) 

p 

being the volume of the parallelepiped spanned by the links issuing from any tetrahe- 
dron vertex. These equalities easily follow from explicit expressions of bivectors in terms 
of link vectors. The consequence is that commutators of such bilinears of the type of 
scalar-scalar and pseudoscalar-pseudoscalar vanish, only those of the scalar-pseudoscalar 
type are nonzero. For example, 

]T F fliy \7r u o vr\ tt° * tt 1 1 = V 2 (F 12 + F 03 - F 31 - F 02 ). (21) 

u>\ J 

Finding A, B, E, F we can find the DB of interest, 

{H»H v } m = {Hi,n^} = {H„H u } + J2^x P [A u ,A,]n x o^ 

\>p 

- £ {H„ ^ * ^}K Xpa ^ - A„)(tt ct o tt*) 

A>p 
a>4> 

+ ^{^,7r A *^}K A ^(A^-A,)(7r CT o7r*) (22) 

A>p 
<T>0 

where we have introduced the notation A^(-) = {H^, •}. The reason for such notation is 
that for the particular case of the solution to the eqs. of motion for Q when R^n^R^ 
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n 



for each \w (i.e. when Q is a true rotation between the local frames) the operator 
coincides with A M (variation at shifting the point \i to the next time leaf). The K^ v \ p 
can be called an inversed to the matrix {tt^ o tt u , 7r A * 7r p }, 



K^ u Xp — Y2V 2 ^ { £ Xpve$p\ + E\pp,£$v £ XpueS p £\ppe$, 



(23) 



(The matrix {tt^ o 71^, 7r A * 7r p } is degenerate, and a strict sense to this term is given by 
solving the eqs. fll8|) , (|19D-) These DB can be further transformed using the following 
two facts. First, note that the operation 



/ - 9 = f ~ E {/> * ^}*W„7r A o tt' 



(24) 



\>P 



on the space of linear combinations of scalars 7r M o -n v projects onto the subspace of such 
combinations of only the squares |vr M | 2 . On these squares = A^ (as consequence of 
the eqs. of motion for (p£). This allows to simplify contribution to {H^,g} which results 
from the dependence of / on the scalars 7r M o tt u . Second, there is some complication 
connected with the (exponential) dependence of on n u @. Suppressing indices, 
consider a function f(ir,Q) of n and Q = fill, II = exp (yj7r + *(p*ir)- The following 
formula can be obtained: 





7T 



7T, un^tu - n^L 
on on 



(25) 



where both derivatives d/dn and Cld/dfl are implied to be antisymmetrised and subscript 
at df/dir denotes the quantity considered as constant at the differentiating. This allows 
to simplify expression for the PB of the two functions f(ir, Oi ) and g(jr, Q 2 ) where Qi = 
Villi, Ui = exp (y^vr + Vi * 7r ) : 



if, 9} = {f,g}n 



71" 



-7T O 



^1 ~„ , ^2' 



niOi^n l5 n 2 ^ 2 ^-n 2 



(26) 



Here {/, p}n are "naive" PB taken in the formal assumption that IT are constants. The 
second term in this formula is responsible for the dependence of IT on n. It vanishes at 
Hi = II 2 as it should because exponential can be eliminated by redefining Q in this case. 

To represent the DB of the Hamiltonian constraints of interest in the form resembling 
that of commutators in the continuum GR (in the HP form) we introduce the notation 



h, 



E 



COS Ot 



(27) 



pi' 



where a pu is angle defect, sin a 



pi' 



np U RiM>/\n>iu>\- The h p arises when we differentiate 



H.^ as complex function over phase variables on which it depends through the angle de- 
fects. Note also that the dependence of Tin on phase variables through the nondynamical 
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ones u v , ip 1 ^, *ip^ need not be taken into account at the calculating the first derivatives 
due to the eqs. of motion for the nondynamical variables. The DB of interest read 



{Hp, Ti u }n{t,) = {hp, h u } a<u + ^ r^jl 71 "* 1 ° ( 



= dhp - dh v 

"A/z^r; < "> 



on 



on 



jrXry 9hp A A dh 



v 



+ 2 [(«fA - tana lyA )A M |n 1/A | - (a M - tana MA )A l/ |n MA |] 

A 

+ i J] ^ X AA» ~ ^A,A /1 )(tt a o vr A ) - £ ^Ap[A^, AJ(vr A o tt") 

Z A A>p 
+ Y, {V ^ * ^}a,n^Ap CT0 [{^, if O 7r*} Qi n - A,(7T CT O TT*)] 

- J2 {hu, ^ * ^,n^w[{ V vr ff o 7^} a , n - A^vr* 7 o tt*)] (28) 

A>p 
cr>(?> 

where IT^ = exp (<fip7r u + V« * 7r ^)> an d the subscript a at the PB means that defect an- 
gles in the denominator in the definition of hp (|27|) are considered formally as constants. 
The first term in the formula obtained is notationally close analog of the PB of the com- 
binations of vector (with coefficients n") and Hamiltonian constraints in the continuum 
GR. 

The DB obtained define then the Jacobian factor entering path integral: 

Det 1 ^ 2 {6,6} = |{Ho, ^l}D( t ){^2, ^3}D(t) + {^0, ^2}D(t){^3, ^i}d(i) 

+ {Ho, W3}d((.){Wi, H2}b(i,)\- (29) 

The rest factors follow by integrating (^-functions in d\i. Integration over d A C, d 6 ip and 
d 6 7i° eliminates 5 A (C) 5 6 (^) S 6 (J2p Integration over d A C reduces S 6 (Z) to 5's of some 
two independent constraints of the type of 7r^ * ix v , fi ^ u. Together with 5 4 ({7r M * 71^ }) 
this gives 8 & {{u a *ti 13 }) (a, f3, j — 1, 2, 3). Finally, the 5 A (ip) annihilate four integrations 
in d 6 ip, and we are left with, e.g., d 2 ip = dip2i dip^ while ipio = ^32, V^o = ^31 = -V'21 - 



^32, 1p. 
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■021- The result reads 



exp li dt 



Y * a o (OA - o ^o) + Y *Px^ x o ^ 



>Det 1/2 {e,e} 



- E ^A m (tt a o tt") 



n 5(n a * n^d 2 ^ w d 6 n a n V 6 flp. 

a>/3 a 11 



(30) 



Some two of the four 5(7Y^)'s can, in principle, be integrated over d 2 ip and the variables 
ippv be excluded as functions of n, Q. Note that despite of that ip seem to enter Hp 
linearly, the implicit dependence on if) through u, if, *(p is much more complex. 
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Due to the Jacobian factor the measure vanishes for symmetrical configurations of 
the system for which the DB {^,W„}d(i) vanish. 

The remaining ^-functions of 7r M * tt u reflect the tetrad structure of bivectors, in 
particular, possibility to place the tetrahedron into 3D space and consider it's area 
vectors instead of bivectors. To pass to the 3D vector formulation in the explicitly 
invariant way one can choose splitting the bivectors and generators of the connections 
into selfdual and antiselfdual parts. To this end we expand 



7T 



ab 



x b +-< b , ^=^Ar/2 (si) 



where ^E^, form basis of (anti-) selfdual matrices obeying algebra of Pauli ones times 
a/— 1. The = { vrf|i = 1,2,3} are vectors in an abstract complex 3D space. The 
generator of the rotation Q = exp uj^ can be decomposed in the same way into + co> At 
and "cu^ which act on the abstract space via vector product — ^uj^ x (•). Denote by ^O^ 
= exp (— ^w^ x (■)) thus obtained representation of Hl^ = exp ( ^lj^) in the abstract 
space. The constraints 7r a * tt 13 = look as 

-jf a . -jfP = +tP ■ +t? /3 . (32) 

This defines ~n a in terms of + n a up to an overall SO(3,C) rotation U so that 

-7? a = Un a , ix a = + 7f a . (33) 

Substituting this into ^ we find that U can be absorbed by ~O m so let us denote 

G„ = -0„U, = +O m . (34) 

The Lagrangian on the constraint surface (kinetic term) takes the form 

L = H \e M ir?(O a O a - OoOo + G a G a - G G ) kl + £ ■ r. (35) 

The Gaussian constraint splits into 

J2(O a -O )n a = 0, (36) 

a 

J2(G a -G )n a = 0. (37) 

a 

Analogously, for curvature matrix R^ u denote by representations of ± R I1U (Rpu = 

+ R^ U ~R liV = ~Rfiv + Rfjiv) in the abstract space. Consider for definiteness i?oi = ^02 
^03- Absorbing U by ~Pqi we find 

P01 == + Poi = exp ( x (.))o 2 3 exp (- + ^tt 3 x (■)), (38) 

Q01 = U-P 01 U = exp(-tpln 2 x(-))G 2 G 3 exp(--tpl7T 3 x(-)), (39) 
% = <±*<- (40) 
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The corresponding contribution to the Hamiltonian reads 



KKi = |n i| arcsin ^l£* (pfeUi + tr Q 01 + Q*Jl + trP 01 ) (41) 

8|n i| v v 

where n m = Uq^ 3 - u^n 2 - tt 2 x 7F 3 . 

Return now to the path integral. Integration of six ^-functions 5(7r a ■ 
over d 3 '^ 1 d 3 ~tx 2 d 3 ~7? 3 results in V~ 2 V 3 U (V 2 = 7F 1 x 7r 2 • 7F 3 ). On the other hand, 
since U is absorbed by "-" components of connection matrices and due to the invariance 
of the measure V 6 tt^ = V 3 + O^V 3 = V 3 O^D 3 G^ the integral over V 3 U decouples 
giving the volume of the gauge subgroup. This leads to the following replacement in the 
measure (|30|): 

n si** * ^) n v^, n ^ =► 1 n d ^ a n ^o,v 3 g, m 

a>/3 ii a a fj, 

The procedure of the derivation of the Ashtekar action by the selfdual-antiselfdual 
splitting of the continuum GR action in the HP form is known in detail (see, e.g., reviews 
|j"5| , |I6|1). We find the following difference of our discrete case from the continuum one: 
now there is no the constraints like the so-called "reality conditions" of the continuum 
GR which would relate and O m . Rather we have a theory with one set of area vectors 
7? M but two sets of connections G^ and O^. The situation is complicated by the fact that 
the dependences of the Hamiltonian on G^ and M does not split, see (f4~l~D. Nevertheless, 
this procedure remains useful, for it allows to remove a part of the gauge degrees of 
freedom in an invariant way. 
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